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Abstract 

Two new QCD sum rules for nucleons in nuclear matter are obtained from 
a mixed correlator of spin-1/2 and spin-3/2 interpolating fields. These new 
sum rules, which are insensitive to the poorly known four-quark condensates, 
provide additional information on the nucleon scalar self-energy. These new 
sum rules are analyzed along with previous spin-1/2 interpolator-based sum 
rules which are also insensitive to the poorly known four-quark condensates. 
The analysis indicates consistency with the expectations of relativistic nu- 
clear phenomenology at nuclear matter saturation density. However, a weaker 
density dependence near saturation is suggested. Using previous estimates 
of in-medium condensate uncertainties, we find M* = 0.64^QQg GeV and 
St, = 0.29to;io GeV at nuclear matter saturation density. 
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Understanding the observed properties of hadrons and nuclei from quantum chromody- 
namics (QCD) is a principal goal of nuclear theorists. The QCD sum-rule approach 
is a particularly useful method for connecting the properties of QCD to observed nuclear 
phenomena ||^-[10|. Recent progress in understanding the origin of the large and canceling 
isoscalar Lorentz-scalar and -vector self-energies for propagating nucleons in nuclear matter 
has been made via the analysis of QCD sum rules generalized to finite nucleon density 
These large self-energies are central to the success of relativistic nuclear phenomenology | ]1 1|| . 

However, the previous sum-rule predictions for the scalar self-energy are sensitive to 
the density dependence of chirally-even dimension-six four-quark condensates. In-medium 
factorization expresses these chirally-even operators in terms of the square of a chirally- 
odd operator whose density dependence is likely to be very different. As such, the density 
dependence of these problematic four-quark operators is generally unknown [0,^0. 

There are various ways to clarify the situation. One approach attempts to better deter- 
mine the density dependence of the four-quark condensates via modeling [O] . One may also 



use other independent information to constrain the in-medium four-quark condensates [jT3 
An alternative approach, which is adopted here, is to derive new QCD sum rules that are 
insensitive to the four-quark condensates. 

In this work, we obtain two new sum rules from a mixed correlator of generalized spin- 
1/2 and spin-3/2 interpolating fields. The spin-1/2 states remain projected, and one gen- 
erates additional sum rules for the nucleon scalar self-energy that are insensitive to the 
four-quark condensates. Our hope is that these new sum rules, along with previous spin-1/2 
interpolator-based sum rules which are insensitive to the problematic four-quark condensates 
I^J^QI, will allow a better determination of the nucleon self-energies in nuclear matter. 

The finite-density QCD sum-rule approach focuses on a correlation function of inter- 
polating fields carrying the quantum numbers of the hadron of interest. The correlation 
function is evaluated in the ground state of nuclear matter instead of the QCD vacuum. 
The appearance of an additional four-vector at finite density, the four-velocity of the nu- 



clear medium, leads to additional invariant functions relative to the vacuum case P-|6|,p!4 



In the rest frame of the medium, the analytic properties of the various invariant functions 
can be studied through Lehman representations in energy. The quasi-nucleon excitations 
(i.e., the quasiparticle excitations with nucleon quantum numbers) are characterized by the 
discontinuities of the invariant functions across the real axis, which are used to identify 
the on-shell self-energies. A representation of the correlation function can be obtained by 
introducing a simple phenomenological Ansatz for these spectral densities. 

On the other hand, the correlation function can be evaluated at large space-like mo- 
menta using an operator product expansion (OPE). This expansion requires knowledge of 
QCD Lagrangian parameters and finite-density quark and gluon matrix elements (conden- 
sates). Finite-density QCD sum rules, which relate the nucleon self-energies in the nuclear 
medium to these QCD inputs, are obtained by equating the two different representations 
using appropriately weighted integrals 0J|,^. 

Consider the correlation function defined by 

nil(g) = 1 1 d'xe^^-{^o\T [xlix)xUo)] l^o) , (1) 
where the ground state of nuclear matter |\l/o) is characterized by the rest-frame nucleon 
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density pn and by the four- velocity u^] it is assumed to be invariant under parity and time 
reversal except for the transformation of u^. The interpolating fields are taken to be |15;16 

Xl{x) = 7^.756'^''= { [ul{x) C75 d,{x)\ u, + P [ulix) C d,{x)] 75 u,{x)} , (2) 



^abc 



{ 



M„ {x)C(jpx db{x 



u^{x)C(7pxUb{x) (T''^7^4(a;)} 



(3) 



where T denotes a transpose in Dirac space, C is the charge conjugation matrix, and /? is a 
parameter allowing for arbitrary mixing of the two independent spin-1/2 interpolating fields 

The correlator of xji and xl gives the sum rules discussed extensively in Refs. In 
particular, the sum rule at the structure, •y^^lu, strongly depends on the even-chirality four- 
quark condensates, {qTiqqTiq)p^ and {qTi XaqqTi Xaq) , where A" are the Gell-Mann 
matrices and Fj are any of the 16 Dirac matrices. The other two sum rules are insensitive 
to these even-chirality four-quark condensates and are given by 



M E,{qq)p, 



+ ^^IJ'^ M'E,{gJia.gq)p,L-''l'' 

7-2/3-5/52 

- Eq {qq)p^ {q^q) 



12 

9 + 10/5 - 29/52 ^_ 



2732 



PJV \"i "i/PJV 
77 



(4) 



and 



Xi St, e 



-(£;?-q2)/Af2 



5 + 2^ + 5/?2 ^,^^ (gtg),^L-^/9 



487r 



+ 



5(5 + 2/5 + 5/?2) _ 

— ^ EqM Eo{qhDoq)p^L ' 



7 + 10/5 + 7/?2 / t r \ r-4/9 
— — ^ M Eo{g^q^a-gq}p^L 



1927r 



+ ^ ^ {^[qUDo^D,q)p, + ^{g.q^a.gq)p,) L"^/^ 



I 5 + 2/5 -|- 5/52 / t \2 r-4/9 
Here Ai denotes the coupling of xi. to the quasi-nucleon state. We have also defined 



(5) 



M: 



N 



Eg = 



q2 + M, 



*2 

N 5 



= - Vq2 + 



r*2 

'AT 5 



(6) 



where and S„ are the scalar and vector self-energies of the nucleon in nuclear matter, 
respectively. The anomalous dimensions of various operators have been taken into account 
through the factor L = \n{M'^ / AI^^b) / Ml^"^ / ^qcd) 0- We have also defined Eo = I - 
^-so/M^ and E'l = 1 — e~*'o/A^^ [sq/M'^ + 1), which account for excited-state contributions 
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In sum rules (^) and (|), contributions proportional to the up and down current quark 
masses have been neglected as they give numerically small contributions. The contributions 
of the even-chirality four-quark condensates to the sum rule (|^) are proportional to current 
quark masses and can thus be neglected safely. Their contributions to (j^) appear only in the 

form of {qioqqioq)pM - j{qY qq7i^q)pN {qio Kqqio Kq)pi, - \{qi^ Kqql^,\aq) p^- 
It is easy to see that these two combinations go to zero in the zero density limit and they 
are proportional to pat in the linear density approximation. Here we will assume in- medium 
factorization for these combinations. To explore the sensitivity of this assumption, we have 
introduced the parameter k in such that deviations from in-medium factorization [k = 
1) may be explored during the Monte-Carlo based uncertainty analysis. The remaining 
nonvanishing four-quark condensates, {qqq'^q)pf^ and {q Xaq q'^ Xaq) p^, have odd chirality 
and vanish at zero density. For simplicity, these condensates are approximated by their 
factorized values. 

The interpolating field xl couples to both spin-1/2 and spin-3/2 states. The three vac- 
uum sum rules obtained from the correlator n^^(g) at ptv = have been studied extensively 
in Refs. p!5|JT6|] . To obtain the finite-density sum rules from this mixed correlator, we follow 
the procedures established in Ref. [0-^. The correlator n^^(g) contains nine distinct struc- 
tures at finite density. Here we only focus on the odd-chirality sum rules at the structures 
7/,7i, and 'jf^^q^, which are 



and 



Here A2 denotes the coupling of xl to the quasi-nucleon state. The contributions of the 
even-chirality four-quark condensates to these two sum rules are multiplied by the current 
quark masses and are thus suppressed. In-medium factorization is adopted for the odd- 
chirality four-quark condensates. Here, dimension-five quark and quark-gluon condensates 
that vanish in vacuum are neglected as their contributions are numerically small [^. The 
other sum rules obtained from the mixed correlator are either dependent on the even-chirality 
four-quark condensates or have vanishing Wilson coefficients for the majority of leading 
order OPE terms. (For example, the Wilson coefficients for the operators, q'j^q, q'jf^D^q, 
■ ■ ■, q'~^^Dy^ ■ ■ ■ D^^^q vanish in this case.) 

In sum rules (^, (^), (|^, and (|), we have also included the dimension-seven condensate, 
{qq)pN (~ ^^)pjv' which was neglected in the previous studies. Inclusion of its contribution 
should make the valid regime in Borel mass larger and the predictions more reliable, as 
shown in the vacuum sum rules ||T6[. As the unfactorized dimension seven operators giving 
rise to {qq)p^ Cl"^) pn ^'^^ chirally odd their density dependence should be qualitatively 



{qq)p. E, L'I^' + 2 E, {qq) p. {q'q)p. L^"' 



1 + 3/? ,_ , .as 



G')p^ L^'^' . 



(7) 



\ {qq)p. E, L'/'^ - %-l {gjla ■ Gq) L"^/^ 



^vr^ - 6471-2 



(8) 
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similar to that of the quark condensate. This behavior arises naturally in the factorized 
form where the density dependence of the gluon condensate is estimated to be a 7% effect. 
There are many other dimension-seven operators; their contributions are assumed to be 
relatively small |]16| and are neglected. 



To analyze the sum rules, we follow the techniques introduced in Ref. |T5[ for deter- 
mining the valid Borel region. We limit the continuum model contributions to 50% of the 
phenomenology, and maintain the contributions of the highest dimensional operators in the 
OPE to less than 10% of the sum of OPE terms. This defines a region in Borel mass M 
where a sum rule should be valid. If there is no region in which both conditions are satisfied, 
the sum rule is considered to be invahd and is discarded. 

A comparison of and (H) indicates that will have a smaller region of validity than 
The leading term of (|^) is proportional to M'* whereas (|^) is proportional to M^. The 

term has greater overlap with excited states and limiting continuum model contributions 
will restrict the upper Borel regime limit. In addition, a comparison of the magnitudes of 
the highest dimension operators in these sum rules indicates the lower limit will be larger 
for (|^ when maintaining the promise of reasonable OPE convergence. It should not be 
surprising to find (J^) to become invalid prior to (P) as density increases. 

QCD sum rules relate the spectral parameters to the condensate values and other pa- 
rameters. Any imprecise knowledge of these QCD inputs will give rise to uncertainties in the 
extracted spectral properties. Here we follow Ref. [0 and estimate these uncertainties via 
a Monte Carlo error analysis. Gaussian distributions for the condensate values and related 
parameters are generated via Monte Carlo. These distributions provide a distribution for 
the OPE and thus uncertainty estimates that are used in the fit. In fitting the sum rules 
taken from the samples of condensate parameters, one learns how these uncertainties are 
mapped into uncertainties in the extracted spectral parameters. 

As in previous works on finite-density sum rules, we use the linear density approximation 
for estimating the in-medium condensates: 

{d)p, = {d),^, + {d)NPN . (9) 

The values of vacuum condensates we use are [|l6l a = —An^ (^Q')o = 0.52 ± 0.05 GeV^, 
b = 47r2 ((a,/7r)G2^o = 1-2 ± 0.6 GeV^ and ml = -{gcW-Gq)o/{qq)o = 0.72 ± 0.08 GeV^. 
The quark mass rriq is chosen to satisfy the Gell-Mann-Oakes-Renner relation, 2mq {qq)o = 



-m: 



fl We adopt cijv = 0.045 ± 0.007 GeV [g, ((a./vr) G^^^ ^ _o,g5 ± g.lO GeV, 
{qHDoq)N = 0.18 ± 0.04 GeV, (^.g^.^g)^ = _o.l ± 0.5 GeV^, {qHDoiDoq)N + 3^(^cgV- 
Gq)N = 0.031 ± 0.010 GeV^ [Q. Since the mixed condensate {gc'qo' ■ Gq)p^ is chirally odd, 
we assume the same density dependence as for the quark condensate. We take /x = 0.5 GeV 
and Aqcd = 150 ±40 MeV and a./vr = 0.0117 ±0.014 at 1 GeV^ |T§. To explore deviations 



from in-medium factorization of the chirally-even four-quark operators contributing to the 
last term of (^) we introduce a 100% standard error and consider k = 1.0 ± 1.0. 

Before proceeding with the numerical analysis it is interesting to examine the density 
dependence of the new sum rules of (J^) and (H). For both of these sum rules, the predomi- 
nant density dependence is governed by the quark condensate. This density dependence is 
common to all terms of these OPEs. As such, the effects of increasing density will be to 
reduce the residue of the pole while the pole position remains largely unchanged. This result 
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should be robust, and is the key result absent in the former in-medium nucleon analysis. The 
approximate invariance of M* + E„ is manifest in (H) and is in accord with the expectations 
of phenomenology. 

Our goal here is to evaluate the degree of consistency between the QCD sum rules and 
the expectations of relativistic mean-field phenomenology. The one firm conclusion from 
previous in-medium studies is that the vector self-energy is positive and a few hundred 
MeV. Hence, we begin by fixing S^, = 0.25 GeV at saturation density and searching for 
a region in which both sides of the QCD sum rules are valid [|1^]. These considerations 
eliminate the sum rules (^) and ([7p from the following analysis as there is no valid Borel 
region as defined above. The implications of this result will be considered in a later work 
20| . Thus we proceed with the sum rules (H) and (|) only. 

The parameter (3 should be selected to minimize continuum model contributions while 
maintaining reasonable higher-dimension operator contributions such that the pole may be 
resolved from the continuum contributions 111^1,0. The sum rule (El) indicates that these 



criteria cannot be satisfied simultaneously. We therefore select /? = —0.7 for this sum rule 
to suppress the continuum contribution from the third term in the OPE. The continuum 
model contribution of the sum rule is independent of (3. We therefore simply set /9 = 0, 
as it is known that the second term of (0) has little overlap with the ground state nucleon 

Ideally, the sum rules of and (1) would be optimized by a six parameter search of the 
two residues, two continuum thresholds and the scalar and vector self-energies. However, 
the final terms of independent of the continuum model are small. A singular value 
decomposition of the covariance matrix reveals that these terms are not identified as a 
degree of freedom. Hence, to proceed, we have no choice but to equate the continuum 
thresholds of (|^) and (^. In practice, this approximation is acceptable when the continuum 
model contributions are restricted as described above ||T6 |. 

Fig. [I| displays the valid Borel regimes for the two sum rules (^) and (^). The corre- 
sponding fit for the central values of the QCD input parameters is illustrated in Fig. |^. 
The distributions for M* and are illustrated in Fig. |. We find M* = 0.70;o:}o GeV 
and = 0.32;{]:?J GeV at nuclear matter saturation density. Normalizing these results to 
the vacuum nucleon mass predicted by (H) provides our final result of M* = 0.64to;o9 GeV 
and = 0.29lHo 

GeV. These results are consistent with the expectations of relativistic 

phenomenology. 

The density dependence of these results near saturation density appears to be weaker 
than previous analyses. Table | summarizes ratios of the finite density spectral parameters 
obtained from (^) and (^) at various densities. Median and standard errors from the median 
are reported as the distributions are not Gaussian. The trend of the ratios is in qualitative 
but not quantitative accord with the expectations of mean-field phenomenology. The un- 
certainties in the ratios for M* + confirm that a determination of binding energies the 
order of 16 MeV is beyond present QCD sum rule analyses. 

The analysis here has focused on two of the four sum rules which are insensitive to the 
even-chirality four-quark condensates and hold promise of providing reliable information on 
the density dependence of spectral parameters. Further examination of the sum rules is 
necessary to determine whether the sum rules which have been found to be invalid are at 
least consistent with the more reliable sum rules. 
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FIG. 1. Illustration of the valid Borel regimes for the sum rules of (^) (large dash) and (fine 
dash). Both continuum model contributions (limited to 50% of the phenomenology) and highest 
dimension operator contributions (limited to 10% of the OPE) are illustrated. 
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FIG. 2. The five parameter fit of the sum rules of (^) (large dash, lower) and (P) (fine dash, 
upper). The QCD— continuum (dashed) curves are hidden by the ground state (solid) curves in 
the near perfect fits. Only 2 of the 51 error bars used in the fit illustrated. 
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FIG. 3. Histogram for S„ and M* obtained from a Monte Carlo sample of 1000 QCD parame- 
ters. 

In summary, we have obtained two new QCD sum rules for nucleons in nuclear matter 
from a mixed correlator of spin- 1/2 and spin-3/2 interpolating fields. These sum rules 
are insensitive to the large and poorly known even-chirality four-quark condensates and 
provide additional information on the nucleon scalar self-energy in nuclear matter. We 
analyzed these new sum rules, along with the sum rules obtained from spin-1/2 correlator. 
Uncertainties of the predictions due to imprecise knowledge of the condensate values and 
related parameters were obtained from a Monte Carlo error analysis. We found that the 
QCD sum-rule predictions for nucleon self-energies are consistent with the magnitudes but 
not the density dependence obtained in relativistic nuclear phenomenology. 



TABLE I. Ratios of the finite density spectral parameters at various densities. Vacuum ratios 
report saturation-density / zero-density results as in S^/Mtv, M*/Mn, (M* + T,y)/M]y, w*/w and 
AIA2/A1A2 while saturation density ratios report finite-density / saturation-density results. 



Self-Energy 


Vacuum Ratio 


Saturation Ratios 




PN = 1.0 


PN = 0.5 


PN = 1.5 




0.30 ± ^oit 


0.88 ± 0.09 


1.09 ± 0.07 


M* 


0.69 ± 


1.06 ±0.05 


0.96 ± 0.05 




0.99 ± 0.03 


1.00 ±0.03 


1.00 ±0.04 


w 


1.01 ±0.04 


1.00 ±0.03 


1.00 ± 0.05 


A? 




0.55 ±0.10 


1.43 ± 0.29 


A1A2 


0.46 ± 0.08 


1.35 ±0.15 


0.72 ± 0.12 
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